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Abstract Systems of coupled oscillators have been seen to exhibit chimera states, 
i.e. states where the system splits into two groups where one group is phase locked 
and the other is phase randomized. In this work, we report the existence of chimera 
states in a system of two interacting populations of sine circle maps. This system 
also exhibits the clustered chimera behavior seen earlier in delay coupled systems. 
Rich spatio-temporal behavior is seen in different regimes of the phase diagram. We 
carry out a detailed analysis of the stability regimes and map out the phase diagram 
using numerical and analytic techniques. 



1 Introduction 

Coupled oscillator systems have long been studied as good models of a variety of 
experimental and natural systems, as well as paradigmatic systems for the obser- 
vation and analysis of complex spatio-temporal behaviour. The study of synchro- 
nised behaviour in systems of identical oscillators, as well as that of oscillators 
with frequencies drawn from a distribution, has resulted in important insights into 
the behaviour of experimental systems, and has also contributed important tech- 
niques for the analysis of extended dynamical systems. In addition to this, a very 
interesting spatio-temporal pattern which is distinct from synchronised behaviour 
was reported by Kuramoto et al. in the case of identical oscillators with symmet- 
rical coupling HI. It was shown that under certain initial conditions and parameter 
values, an array of identical oscillators split into two groups: one showing coher- 
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Fig. 1 Snapshots of the oscillators of the two populations in (a) two clustered state (b) clustered 
state (c) clustered chimera (d) chimera. The blue lines (thin line) correspond to group 1 where all 
the oscillators are synchronised. In (c) the oscillators from group 2 overlap with group 1 . The order 
parameter for (e) the completely synchronized state (f) aperiodic chimera (g) stable chimera. 



ent behavior and another incoherent behavior. Such states were named as chimera 
states by Strogatz et. al 13. The splitting of globally coupled oscillator population 
with distributed frequencies into synchronized and desynchronized parts was earlier 
reported |3]|4]|5|- Various theoretical and numerical results have been reported since 
this discovery J6] 13 . An exact result for the stability and bifurcations of a system of 
two interacting populations of oscillators has been obtained J8). It was shown that 
in a globally coupled network of oscillators with delay feedback, the single spatially 
connected region was replaced by a number of spatially disconnected regions of co- 
herence with regions of incoherence in between. These states were named clustered 
chimera states(|9] [10). Chimera states have been described as the natural link be- 
tween coherent and incoherent states ifTTl . The spatiotemporal pattern for the coher- 
ent and incoherent regions forming the chimera states has been studied Ifl2l . In this 
work, we study the existence and stability of chimera states in two interacting pop- 
ulations of coupled sine circle maps. These systems are the map analogs of coupled 
oscillator systems. Coupled map systems show many of the varied phenomena ob- 
served in the case of continuous extended systems, but are far more numerically and 
analytically tractable. Coupled map systems have been seen to show synchronised 
solutions, spatiotemporal intermittency, spatiotemporal chaos |[T3l[T4l . However, to 
the best of our knowledge, chimera states have not been reported earlier in coupled 
map lattice systems. 
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2 The model 



The single sine circle map evolves via the evolution equation 

e,+i =f(e t ) = e t +Q--^-sm(27te t ) modi (i) 

Here, K is the nonlinearity parameter and Q is the winding number of a single 
sine circle map in the absence of the nonlinearity. This map shows a tendency to 
mode lock as the parameter K is increased and the phenomena of Arnold tongues 
are seen in the K — Q space lfT31 . The winding numbers of the mode locked tongues 
show a Devil's staircase structure. There are several studies on the spatiotemporal 
dynamics of diffusively coupled sine circle map lattices on regular sites [161 . The 
model that we consider consists of two interacting populations of sine circle maps 
with every element of one population coupled via a coupling parameter £j to all the 
other elements in the same group and via another parameter £2 to all those in the 
other group. The evolution equation is given by 

C 1 (/) = (2- ei - e2 )(e;(/) + r2-^sin(2^(/))+ £ fspnf^) 

<r'=l ls o> j= i 

+ Q. sin(27r0„ CT '(y')) modi (2) 

2n 

where £n = £22 = £1 and £12 = £21 = £2- 0(0 is the angle at time t and lies between 
and 1 . The parameters Q and K are taken to be uniform at each site. The values of 
£1 and £2 should lie between and 1 and £1 + £2 = 1. We considered 1024 oscillators 
in each group and the results remain the same with fewer oscillators. The first popu- 
lation is given identical initial conditions while the second set is given random initial 
conditions between and 1 . The system exhibits rich spatial dynamics depending 
upon the value of K and £2 for this initial condition. The most common solution 
for large values of K are the completely synchronized state and the two clustered 
state, where all the oscillators of each group have the same value of the phase but 
the phase values for each group differ from each other are shown in Fig. 02 a )- As 
the value of K decreases the oscillators of group two bifurcate to the clustered state 
as shown in Fig[TJb). With further decrease in the value of K all the oscillators of 
group two completely desynchronise, which correspond to the chimera state seen in 
FigDId). One interesting solution that we observed was the clustered chimera state 
where the group two oscillators again separate into a completely synchronized state 
while the rest of the population remains desynchronized Figfllc). In Fig. [Tie) half 
of the population of the second group has the same state value as the first group 
and hence overlaps and the rest remains distributed over the phase space. For all the 
values of K and Q the oscillators of group one remain completely synchronized. 

The synchronization within the population of each kind can be characterised by 
the order parameter r(t) = |(exp(/0y(f)))2| where the angle brackets denotes the 
average over the oscillators having the same dynamics. Thus for the completely 
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Fig. 2 (Color online)The 
parameter space plot for 
£i = 0.9. The red (■ )-s show 
chimera states, black (+)- 
s show clustered chimeras, 
cyan (A)-s show clustered 
states, blue (o)-s show three 
clustered states, black (*)-s 
show two clustered states and 
yellow (o)-s show complete 
synchronous states 




synchronized, and the chimera state, the average is taken over all the oscillators of 
group 2. The order parameter r will be 1 for a completely synchronized state as the 
phases of all the oscillators are the same in this state, and will take values between 
and 1 when there is desynchronization or clustering. When r takes a constant value 
the corresponding state is a stable state. A time varying r corresponds to either a 
breather solution or an unstable state. Chimera states can be statistically stable |fl~) 
in some cases leading to stable order parameters, in other cases the order parameter 
for chimera states can be oscillating or aperiodic (8). In our case, we had stable 
chimera states which bifurcate to the clustered state through an unstable aperiodic 
chimera state. Fig. [TJe) shows the order parameter for the group 2 oscillators for 
the completely synchronized state and FigsQlf) an d (g) shows the aperiodic and 
stable chimera respectively. For the clustered chimera and the clustered state there 
is further division in the population and r has to be defined separately for each group. 

We explore the phase diagrams of the system at K — £2 — e values which exhibit 
the variety of possible solutions and also provide clues to the bifurcation sequences 
that take place. 



3 The phase diagrams and stability analysis 

3.1 The K — fl phase space 

We explore the K — Q phase space with £\ = 0.9, where a good spectrum of dy- 
namical behaviour is seen and map out the observed types of spatial behavior. The 
value of K is varied from to 1 and Q from to 0.5. The phase diagram is sym- 
metric about £2 = 0.5, hence values beyond Q = 0.5 have not been plotted. Figf2] 
is the K — Q parameter space plotted for £\ = 0.9. Low values of K favor chimera 
states whereas for high values of K the system exhibits ordered spatial behavior. For 
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values of K near 1 the system either is either in a completely synchronized state or a 
two clustered state depending upon the Q values for the initial conditions mentioned 
before. The yellow (o)-s and the black (*)-s represent the completely synchronised 
and two clustered state respectively and form the most observed states. Thus for a 
majority of parameter values the system is spatially ordered and forms a symmet- 
ric pattern about Q = 0.25. As the value of K is decreased for values of Q near 
and 0.5 the oscillators of group 2 bifurcate into a three clustered state —5- clustered 
chimera — > chimera state. In the remaining regions the clustered state bifurcates to 
a chimera state. Here the order parameter of the chimera state is aperiodic near the 
clustered state and then forms a stable chimera. The temporal behavior throughout 
the Q —K space is chaotic except for the completely synchronized region. 

Now we study three e — £2 sections at three distinct values of K. One is the shift 
map case where K = 0, another is the K = 1 plane which is known to mode lock and 
the third is a low value (K — 0.1), where chimera states are seen and there is rich 
dynamical behavior. 



3.2 The Shift Map case 

The simplest case of sine circle maps is the shift map case with K = 0. For a single 
shift map given by 

9,+ i = 9, + Q modi (3) 

the system has periodic orbits for rational values of Q . The corresponding equation 
for the coupled system is 

0»+i(O = (tf(O + fl)+ t TrXX'C0 + <2) modi (4) 

a'=\ j=\ 

For the shift map case, the system always evolves to a chimera state for the initial 
conditions mentioned earlier. The numerically evaluated largest lyapunov exponent 
shows that the system is chaotic in the entire £2 — E\ phase space. The order parame- 
ter defined by r(t) — |(exp(/0 ; (f))}2| remained constant with very small fluctuations 
which were due to the finite size effects. We checked this by varying the number of 
nodes in the constituting groups. Thus the shift map case has spatially stable chimera 
states that are temporally chaotic. The Jacobian matrix for the coupled system may 
be written as 




where 
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Fig. 3 (Color online)A) The parameter space for the shift map case for the coupled sine-circle 
maps. Chimera states exists for all values of Q&£[ . B)The parameter space plot for K = 1 . The 
red (A)-s show completely synchronized fixed point solutions, black (*)-s show completely syn- 
chronized solutions, blue (°)-s shows two clustered chaotic solutions. For K = 1, there are no 
chimeras. 



/(2-^ £l -e 2 )f'(9,](l)) 



A = 



1/(^(2)) 



e if'(e,}(N)) 



E if'(oKN)) \ 



(2-^-^/(^(2))/ 



D 



(2-^i ei - £2 )^(0„ 2 (l)) 



fg'(0 2 (l)) (2-^l ei -e 2 y(0„ 2 (2)) 



&?'(#(#)) 



(2-Vei-£2V(0„ 2 (2)). 



5 



C 



/tm 2 (i)) i^'(0„ 2 (2)) 

^'(0 2 (l)) |^(0„ 2 (2)) 
VtW, 2 (1)) l^'(0, 2 (2)) 
/ft/W)) |/'(^(2)) 



Vf/'(^(1)) |/'(0„ 1 (2))...|/'(9» 1 (iV))/ 



%'(0„ 2 (A^))\ 
|s'(0„ 2 (AO) 



e is'{e 2 n (N))J 
|/'(0 n 1 (^v))\ 



We have f'(6*(N)) = 1 - Kcos(27t9^(N)) and g'(0„ 2 (AO) = 1 - ^cos(27T0 2 (iV)). 
The eigen values need to be numerically obtained since the stable state is the 
chimera state. 
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Fig. 4 A ) The largest Lyapunov exponent obtained numerically for K = 1.0. B) LE evaluated 
from Eq. [5] for K=\. The blue (o)-s represent the periodic solutions and cyan (A)-s are the 
chaotic solutions 



3.3 Linear Stability Analysis for K —I 



FigfS] shows that in the £i — £2 phase space plane with K = 1, the system is always 
having either a complete synchronized state or a two clustered state. We study this 
plane analytically to obtain the bifurcation boundary. In the completely synchro- 
nized state 0„' (i) = 6%(i) = 0* where i varies from 1 to N. In this case the Jacobian 
matrix reduces to 
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This can be reduced by a similarity transformation to the form 
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J'=.f'(0n) 



A'+B' 
A'-B 1 



Now the Jacobian matrix has a block diagonalized form with each nonzero block 
being a circulant matrix such that 
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This matrix has eigenvalues 2f'(9*), 2e\f'(d*) and (N -2)- fold degenerate 
eigen values f(9*). The Lyapunov exponents in terms of the eigen- values of the 
Jacobian matrix can be written as 



A, = Ilim£ln|2/'(e;)| 
A 2 = ilim£ ln|2ei/'(e;) 

T T— >oo 



t=\ 



and 



A 3 = A4 = ... = X N = i lim £ln | 



(5) 
(6) 

(7) 



The numerically obtained value of the largest lyapunov exponent for the syn- 
chronised state agrees with the value given by EqO FigHJa) shows the numerical 
LE for the entire £2 — £1 and the system is chaotic for most of the values. In the peri- 
odic regions, for £2=0 and 0.5 the system has fixed point solutions and is periodic 
otherwise as seen in the figure. 



3.4 The Q.—£\ parameter space for K = 0.l 

The third slice along K = 0.1 exhibits rich spatial dynamics. As noted earlier, 
chimera states are easily found at the lower values of K. For 6] above 0.5 and for £2 
near 0&0.5 the system has clustered chimera states (black +-s) which bifurcate to 
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Fig. 5 (Color online) A) £2 — E] plot with K = 0.1. Red (-)-s are the chimeras and black (+)-s 
are the clustered chimeras.The cyan (A)-s shows clustered states, blue (o)-s shows three clustered 
states, black (*)-s shows two clustered states and yellow ( A)-s shows complete synchronous states. 
B) The largest Lyapunov exponent shows that the system is always chaotic. 



clustered states (cyan A) as shown in Figj5] . For fairly large regions in the param- 
eter space the system exhibits chimeras (red •). The order parameter corresponding 
to the chimeras are by and large constants and hence stable chimeras are seen in this 
regime. The system also has two clustered solutions in this region which bifurcate 
to the unstable chimeras and then to stable chimeras. For values of £i less than 0.5 
in addition to these solutions the system also has a completely synchronized region 
(yellow A). Thus, we note that increasing the nonlinearity parameter K tends to re- 
duce the region which is available for the chimera solutions. A 3 — d phase diagram 
will illustrate this aspect further. We plan to explore this in future work. 



4 Conclusions 

In this paper, we proposed a couple map model with two interacting species of sine 
circle maps coupling. Chimera states are easily seen in this system given two sets 
of initial conditions for the two species, one random and one synchronised. Special 
choices of initial conditions are not required. The K = i.e. the shift map version of 
this models shows temporally chaotic chimeras all over the parameter space. The in- 
troduction of K stabilises non-chimera states , but chimera states are still easily seen 
at low values. Clustered chimera states where the desynchronized population again 
splits into a synchronized and desynchronized population are seen here, without the 
introduction of delays. The stability analysis for synchronised states is carried out 
analytically and numerically, so that bifurcations from the synchronised values can 
be obtained. We observe that the stable chimera bifurcates to the clustered state 
through an aperiodic chimera. The full 3— dimensional phase diagram remains to 
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be explored, hence other bifurcation sequences for the chimeras may be possible. 
We hope to explore these directions in future work. 
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